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ABSTRACT 


An  investigation  of  natural  convection  heat  transfer  from 
a  flush  mounted  heater  on  a  larger  horizontal  substrate  in 
water  has  been  conducted.  The  focus  of  the  present 
investigation  was  on  the  buoyant  flow  above  the  heat  source. 
Temperature  measurements  in  the  buoyant  plume  were  made  at 
various  locations  and  heights,  for  a  number  of  input  powers 
using  a  traversing  thermocouple  probe.  Simple  signal 
processing  was  performed  on  the  temperature  measurements  to 
attempt  to  investigate  the  reported  meandering  motions  in  the 
plume  and  the  transition  from  laminar  to  turbulent  buoyant 
flow.  These  measurements  are  interpreted  based  on  flow 
visualizations . 
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I .  INTRODUCTION 

A.  NATURAL  CONVECTION  ABOVE  A  HORIZONTAL  HEATED  FLAT  PLATE 

Natural  convection  heat  transfer  from  flat  surfaces 
continues  to  be  an  area  of  study  that  requires  further 
research.  Many  scientific  and  engineering  applications,  such 
as  meteorological  thermal  rises  or  plumes  emanating  from  the 
Earth's  surface,  are  driven  by  natural  convection  heat 
transfer.  Investigations  of  this  mode  of  heat  transfer  from 
flat  surfaces  have  been  conducted  for  a  variety  of  planforms 
and  orientations,  and  using  a  variety  of  techniques. 

Husar  and  Sparrow  [Ref.  1]  researched  the  nature  of  free 
convection  flow  fields,  emphasizing  the  influence  of  the 
heated  area  on  the  shape  of  the  fields  near  the  surface.  They 
showed  that  the  flow  is  normal  from  of  the  edge  of  the  surface 
towards  the  center.  Depending  on  the  planform,  bisectors  are 
created  by  the  colliding  flow  paths  as  the  boundary  layer 
develops  into  an  ascending  plume.  Ackroyd  [Ref.  2]  developed 
a  numerical  analysis  based  on,  and  in  good  agreement  with,  the 
experimental  results  obtained  by  Husar  and  Sparrow. 

Goldstein  et.  al.  [Ref.  3]  investigated  the  natural 
convection  mass  transfer  adjacent  to  a  horizontal  downward¬ 
facing  plate.  Using  a  naphthalene  sublimation  technique,  a 
common  correlation  for  horizontal  plates  with  specific 
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characteristic  lengths  in  the  Sherwood  and  Rayleigh  numbers 
was  obtained.  The  results  were  compared  to  a  number  of  other 
studies  and  some  differences  were  noted  due  to  deviations  in 
geometry  and  experimental  technique.  Lloyd  and  Moran  [Ref.  4] 
used  an  electrochemical  technique  to  measure  mass  transfer. 
They  developed  correlations  for  turbulent  and  laminar  mass 
transfer,  the  latter  agreeing  with  Goldstein  et.  al. 

In  the  field  of  meteorology,  natural  convection  in  a  large 
body  of  fluid  leading  to  freely  rising  plumes  is  of  great 
concern.  This  led  Pera  and  Gebhart  [Ref.  5]  to  conduct 
analysis  and  experiment  for  horizontal  and  slightly  inclined 
surfaces.  Using  interferograms,  a  comparison  was  made  between 
analytical  and  experimental  results  which  were  found  to  be  in 
fair  agreement  depending  on  Prantl  and  Grashof  numbers.  Pera 
and  Gebhart  [Ref.  6]  then  repeated  their  investigation  to 
include  small  disturbances.  Again,  using  interferograms, 
their  goal  was  to  compare  an  analysis  of  the  stability  of  flow 
subjected  to  small  disturbances,  to  experimental  results  in 
order  to  provide  greater  insight  into  instability  and  flow 
transition. 

Al-Arabi  and  El-Reidy  [Ref.  7]  studied  natural  convection 
over  heated  isothermal  plates.  Investigating  a  number  of 
planforros,  they  developed  new  local  and  average  heat  transfer 
correlations  for  both  laminar  and  turbulent  boundary  layers. 
Compared  to  equations  developed  by  Fishenden  and  Saunders,  the 
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results  obtained  for  the  laminar  case  were  30%  higher  and  for 
the  turbulent  case  11%  higher. 

Natural  convection  from  a  horizontal  flat  heater,  flush 
mounted  in  an  unheated  substrate,  was  studied  by  Hickey  [Ref. 
8].  Hickey  investigated  the  effects  of  both  steady  and 
pulsatile  input  powers  on  the  surface.  A  correlation  for  the 
resulting  natural  convection  heat  transfer  was  developed. 
From  the  data  obtained  for  pulsatile  input  powers  it  was  noted 
that  the  effect  of  pulsations  resulted  in  minimal  heat 
transfer  enhancement  as  compared  to  similar  steady  power 
inputs . 

B.  BUOYANT  PLUME  ABOVE  HEATED  HORIZONTAL  SURFACE 

Free  boundary  flows  such  as  buoyant  plumes  or  jets, 
emanate  from  heated  surfaces  due  to  density  gradients 
generated  by  the  differences  in  temperature.  The  density 
gradient  produces  proportional  body  forces  in  a  gravitational 
field  inducing  free  convection  currents.  These  flows  have 
been  the  subject  of  continuing  research. 

Spalding  and  Cruddace  [Ref.  9]  performed  a  theoretical 
analysis  for  the  steady  laminar  buoyant  flow  above  a  line  heat 
source.  Their  analysis  was  followed  by  Fujii  [Ref.  10]  who 
also  studied  the  flow  above  a  point  heat  source.  Forstrom  and 
Sparrow  [Ref.  11]  applied  these  analyses  to  experiments  they 
performed  above  a  heated  wire.  They  observed  plume 
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oscillations  and  determined  the  transition  and  onset  of 
turbulence  in  terms  of  a  modified  Grashof  number. 

A  study  of  the  transition  of  plane  plumes  was  conducted  by 
Bill  and  Gebhart  [Ref.  12].  Their  investigation  in  air 
clarified  some  of  the  mechanisms  involved  in  the  transition  to 
turbulent  buoyant  flow.  Evaluating  not  only  quantitative 
data,  but  also  visual  criteria  from  interferograms,  they 
identified  the  effects  of  disturbances  in  the  transition 
process  and  were  able  to  define  the  limits  of  the  transition 
regime  in  air. 

C.  PRESENT  STUDY 

The  present  study  continued  the  investigation  of  natural 
convection  heat  transfer  from  a  flush  mounted  heater  on  a 
larger  horizontal  substrate  in  water  begun  by  Hickey  [Ref.  8], 
The  primary  focus  was  on  the  buoyant  plume  emanating  from  the 
surface.  While  previous  studies  involved  point,  line,  or 
horizontal  finite  sources,  little  is  known  about  free  boundary 
flow  from  a  finite  source  surrounded  by  an  unheated  area  such 
as  a  substrate  mounted  electric  component. 

The  objectives  of  this  study  were: 

•  To  fabricate  a  traversing  thermocouple  probe  and  to  use 
the  probe  to  measure  the  temperature  distribution  within 
the  buoyant  plume. 

•  To  investigate  meandering  motions  of  the  plume  reported  in 
earlier  studies. 

•  To  investigate  the  transition  of  the  free  boundary  flow 
above  the  heat  source. 
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•  To  visualize  the  buoyant  flow  above  the  surface  in  order 
to  interpret  the  meandering  and  transition  process  of  the 
flow. 
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ZZ.  EXPERIMENTAL  APPARATUS  AND  PROCEDURES 


A.  EXPERIMENTAL  APPARATUS 

The  experimental  apparatus  consisted  of  four  primary 
interrelated  systems  and  two  independent  secondary  systems  as 
seen  in  Figure  1.  The  primary  systems  included:  the  power 
supply  system;  the  test  surface  heater  system;  the 
thermocouple  probe  traversing  system;  and,  the  data 
acquisition  and  reduction  system.  The  secondary  systems 
included:  the  flow  visualization  system  and  the  water 
deaeration  and  filtering  system.  The  power  supplied  to  the 
heater,  the  temperature  measurements  on  the  test  surface  and 
in  the  buoyant  plume,  the  position  of  the  thermocouple  probe, 
and  the  temperature  of  the  ambient  fluid,  are  the  related 
variables  of  the  primary  systems  that  the  experiment  depended 
upon.  Flow  visualization  and  water  deaeration/purification 
were  independent  systems. 

All  of  the  systems  mentioned  above  are  described  by  Gaiser 
[Ref.  13],  Haukenes  [Ref.  14],  Akdeniz  [Ref.  15],  Larsen  [Ref. 
16],  and  Hickey  [Ref.  8];  described  are  the  basic  set-up,  the 
water  tank,  data  acquisition  system,  and  water  filtration 
system,  and  additions  or  modifications  such  as  the  deaeration 
system  or  the  thermocouple  traverse  system.  The  following 
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input  power 


summarizes  those  systems  and  describes  modifications  which 
reflect  the  current  condition  of  the  apparatus. 

1.  Test  Surface 

The  test  surface  was  described  in  detail  by  Hickey 
[Ref.  8].  The  surface  consisted  of  a  30.48  cm  x  30.48  cm 
horizontal  substrate  with  a  9  cm  x  9  cm  heater  flush  mounted 
in  the  center.  Embedded  in  the  surface ,  under  the  heater, 
were  twenty-five  0.013  cm  copper-constantan  thermocouples  as 
seen  in  Figure  2.  The  heater,  which  provided  a  uniform  heat 


1.5 

2.0 

2.0 

2.0 

2.0 

1.5 


Figure  2.  Thermocouple  Locations  under  Test  Heater. 


flux,  was  a  gold/tin  oxide  film  with  a  surface  resistance  of 
2.6  ohms/square;  net  thickness  of  film  and  substrate  was  0.165 
mm.  The  film  was  attached  above  the  thermocouples  using  a 
thermally  conductive  epoxy.  The  power  leads  were  attached  to 
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the  film  using  a  silver  loaded  "ink"  with  a  resistance  of  0.2 
ohms/square.  When  the  power  leads  were  attached,  the  entire 
surface  had  a  resistance  of  approximately  3  ohms.  Covering 
the  entire  test  surface  was  a  thermochromic  liquid  crystal 
(TLC)  sheet.  Used  to  ensure  a  smooth  flat  surface,  the  sheet 
also  provided  a  visualization  of  the  temperature  patterns  on 
the  surface;  it  was  calibrated  in  water  and  showed:  yellow  at 
approximately  41°C,  green  at  43°C,  blue-green  at  45°C,  and 
dark  blue  at  temperatures  above  47 °C. 

Because  the  thermal  properties  of  the  various 
components  of  the  test  surface  were  different,  care  had  to  be 
taken  when  increasing  the  power  to  the  heater  because  sudden 
thermal  expansion  would  cause  it  to  fail.  This  occurred  when 
the  gold/tin  oxide  film  heater  was  rapidly  heated  and 
conduction  to  the  substrate  was  slow;  the  gold/tin  oxide  would 
separate  from  the  film  upon  which  it  was  attached  because  of 
the  bonding  power  of  the  epoxy.  Reattachment  of  a  new 
gold/tin  oxide  heater  proved  difficult  after  failure  occurred 
and  required  slight  modification  of  the  surface. 

The  difficulty  in  rebuilding  the  test  surface  was  with 
the  attachment  of  the  power  leads  to  the  new  gold/tin  oxide 
film.  The  original  surface  was  built  from  the  top  down  and 
attachment  of  the  power  leads  to  the  gold/tin  oxide  film  prior 
to  affixing  the  film  to  the  surface  was  possible.  In  the 
latter  case,  the  power  leads  could  only  be  affixed 
simultaneously  with  attachment  of  the  film  to  the  substrate. 
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Because  the  loose  wires  of  the  power  leads  were  difficult  to 
work  with,  copper  strips  were  soldered  to  the  power  leads,  and 
the  strips  embedded  in  the  substrate  as  seen  in  Figure  3 . 


Figure  3.  Bottom  View  of  Test  Beater. 

With  a  solid  conductive  strip,  the  gold/tin  oxide  film  could 
easily  be  attached  to  the  surface  by  spreading  the  epoxy  on 
the  substrate,  the  silver  loaded  "ink"  on  the  copper  strips, 
and  pressing  the  film  in  place.  The  result  was  that  a  3  ohm 
surface  resistance  was  easily  obtained  across  the  rebuilt 
heater. 
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2.  Power  Supply  System 

The  entire  experiment  is  driven  by  the  power  supply 
system  which  provided  the  power  to  the  test  surface  heater. 
The  driver  in  this  system  was  a  HP-85  micro  computer  which  was 
programmed  to  control  various  power  inputs  to  the  test  surface 
heater.  The  computer  directed  a  HP-59501B  Isolated  DAC/Power 
Supply  Programmer,  which  in  turn  controlled  the  HP-6286A  D.C. 
Power  Supply  that  actually  provide  the  power  to  the  heater. 
As  shown  in  Figure  4,  the  heater  was  connected  in  series  with 


HP-59501  B  HP-6286A  D.C. 
Programmer  Power  Supply 


Figure  4.  Power  Supply  System. 

a  0.101  Ohm  precision  shunt.  By  obtaining  the  voltage  across 
the  precision  shunt  (V*,,*)  and  across  the  heater  ,  and 
by  knowing  the  resistance  through  the  precision  shunt  , 
the  input  power  to  the  surface  can  be  computed  using  the 
following  equation: 
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Powez  =  (-^2^)  VMC„ 
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3.  Thermocouple  Probe  Traversing  System 

The  thermocouple  probe  traversing  system  was  used  to 
obtain  temperature  measurements  at  specific  locations  in  the 
buoyant  plume.  The  system  consisted  of  a  thermocouple  probe, 
a  traversing  system  with  three  degrees  of  freedom,  and  a  IBM 
PC  clone  "notebook"  computer. 

The  thermocouple  probe  was  made  by  passing  the  ends  of 
0.0076  mm  copper-constantan  thermocouple  through  two  bent 
glass  capillaries  which  were  bonded  together.  The  ends  of  the 
0.0076  mm  thermocouple  wires  were  then  welded  to  0.76  mm  wire. 
The  0.76  mm  wires  and  glass  capillaries  were  then  pulled 
through  a  stainless  steel  tube,  bonded  in  place,  and  the  ends 
of  all  the  tubes  were  sealed.  The  probe,  shown  in  Figure  5, 
was  then  attached  to  the  traverse  via  a  stainless  steel  rod. 

The  traversing  slides,  manufactured  by  VELMEX, 
consisted  of  two  motor  driven  slides  and  one  manually 
positioned  slide.  The  powered  slides  were  positioned  by  a 
86MM-2  two  axis  stepping  motor  controller/dr iver.  These  were 
used  to  move  the  probe  within  a  horizontal  plane  according  to 
position  commands  entered  into  the  "notebook"  computer.  The 
manual  slide,  with  a  30  cm  vernier  scale,  was  used  to  move  the 
probe  in  specific  planes  above  the  heater  surface.  The  entire 
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clone  personal  computer.  The  HP  computer  was  programmed  to 
interrogate  the  data  acquisition  unit  for  monitored 
temperature  and  voltage  measurements ,  and  to  transmit  the 
reduced  information  to  the  IBM  clone  via  an  export  program. 
The  IBM  clone  was  used  to  process  the  imported  data  through 
spreadsheet  and  graphics  programs. 

S.  Secondary  Systems 

Not  directly  linked  to  the  four  primary  systems,  the 
flow  visualization  and  deaeration/ filtration  systems  were 
considered  secondary.  The  flow  visualization  system  consisted 
simply  of  a  He-Ne  laser  whose  beam  was  spread  into  a  sheet 
through  a  cylindrical  lens  to  illuminate  suspended  particles 
in  the  water.  As  seen  in  Figure  6,  a  SLR  camera  was  used  to 


Camera 


J  1  Laser 


Test  Surface 


Cylindrical  Lens 


Figure  6.  Flow  visualization  System. 
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photograph  the  flow  path  of  the  particles  through  the  buoyant 
plume  as  they  passed  along  the  sheet  of  generated  laser  light. 

The  deaerating  unit  consisted  of  a  deaerating  tank,  a 
vacuum  pump,  a  heater,  and  a  discharge  pump.  With  the  vacuum 
pump  lowering  the  pressure  in  the  deaerating  tank  to  68  to  85 
kPa,  water  was  suctioned  into  the  tank  through  spray  nozzles. 
With  the  pressure  low  and  the  heater  increasing  water 
temperature,  entrained  air  escapes  from  the  water  as  it 
settles  into  the  tank.  The  deaerated  water  is  then  returned 
to  the  test  tank  via  the  discharge  pump.  Attached  to  the 
deaerating  unit  return  line  is  the  filtration  unit.  This  unit 
suctions  water  from  the  test  tank  and  passes  it  through  four 
cartridge  filters  which  remove  particulates  and  deionize  the 
water. [Ref.  13] 

B.  EXPERIMENTAL  PROCEDURE 

The  experimental  program  consisted  of  three  phases:  test 
surface  temperature  evaluation,  buoyant  plume  temperature 
evaluation,  and  visualization  of  the  buoyant  flow.  In 
preparing  for  all  phases,  a  mechanical  mixer  was  operated  for 
at  least  five  minutes  to  dissipate  temperature  stratification 
within  the  test  tank.  All  primary  systems  were  then  activated 
while  the  tank  was  permitted  to  reach  quiescence,  which  took 
approximately  30  minutes. 


15 


l.  Test  Surface  Temperature  Evaluation 


The  procedure  for  evaluation  of  the  test  surface  was 
completely  described  by  Hickey  [Ref.  8],  Power  was  provided 
to  the  test  heater.  After  approximately  15  minutes,  when 
steady  state  was  reached,  temperature  measurements  were  taken 
on  the  surface  of  the  heater  over  a  nine  minute  period.  This 
was  done  for  seven  different  constant  and  two  different 
periodic  power  levels. 

2.  Buoyant  Plume  Temperature  Measurements 

Following  the  surface  temperature  measurements  for  a 
specific  power  level,  temperature  measurements  in  the  buoyant 
plume  were  obtained.  Measurements  were  obtained  within  one 
quadrant  over  the  heater  surface  at  nine  locations  and  six 
elevations. 

To  obtain  the  plume  temperature  measurements,  the 
thermocouple  probe  was  used.  Measurements  always  commenced 
with  the  probe  in  the  center  position  (45  mm,  45  mm)  at  an 
elevation  of  2.0  mm  above  the  surface.  At  that  position,  the 
data  acquisition  program  was  run  and  1300  temperature 
measurements  were  obtained  over  a  period  of  260  seconds. 
Early  evaluation  of  the  data  acquisition  program  and 
measurements  showed  that  a  5  He  sampling  frequency  was 
adequate  to  meet  Nyquist  sampling  criterion. 

When  plume  temperature  measurements  at  one  position 
were  obtained,  the  probe  was  repositioned.  This  was 
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accomplished  by  typing  movement  instructions  on  the  “notebook" 
computer  which  relayed  these  instructions  to  the 
controller/driver  and  caused  the  traverse  slides  to  move  to 
the  new  position.  After  the  plume  was  given  time  to  settle 
following  disturbance  by  probe  positioning,  the  above 
procedures  were  repeated  until  readings  were  obtained  for  the 
nine  positions  at  six  different  elevations. 

3.  Flow  Visualisation 

The  final  phase  of  the  experiment  program  involved 
flow  visualization.  With  all  temperature  measurements 
obtained,  contamination  of  the  test  tank  with  the  flow 
visualization  particles  was  possible. 

With  the  test  heater  powered  at  specific  levels, 
ground  pliolite  powder  was  dispersed  in  the  test  tank.  These 
particles  are  easily  suspended  in  the  water  (specific  gravity 
1.026)  and  scatter  laser  light  when  illuminated.  A  vertical 
plane  of  laser  generated  light  was  passed  over  the  center  of 
the  test  heater  to  illuminate  the  particles  which  could  then 
be  photographed  at  extended  time  exposures.  Photographs  were 
taken  for  15,  30,  60,  and  120  second  exposure  times  at  a  low, 
medium,  and  high  input  power.  The  procedure  was  described  in 
detail  by  Gaiser  (Ref.  13]. 
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III.  ANALYSIS  OF  FLOW  VISUALIZATION 
Flow  visualization  was  conducted  for  a  number  of  input 
powers.  The  objective  was  to  obtain  representative  flow 
patterns  for  a  low,  medium,  and  high  heater  power  level. 
Results  for  five  different  power  levels,  of  which  one  was  also 
the  average  of  a  pulsatile  triangular  wave  input  power,  were 
successfully  obtained.  Photographs  were  taken  of  two 
perpendicular  central  planes  above  the  test  surface. 

The  focus  of  the  flow  visualization  program  was  to  verify 
some  of  the  reports  made  in  earlier  studies.  Forstrom  and 
Sparrow  [Ref.  11]  observed  that  for  purely  laminar  flow  the 
plume  exhibited  slow,  regular  meandering  motion.  They  also 
noted  that  depending  on  the  heating  rate  and  elevation  above 
the  source  the  onset  of  transition  was  indicated  by  turbulent 
bursts.  Bill  and  Gebhart  [Ref.  12]  observed  the  effects  of 
disturbances  in  aiding  the  onset  of  transition. 

All  the  photographs  were  taken  with  various  exposure  times 
in  the  range  15  to  120  seconds.  Because  of  that,  the  flow 
paths  of  the  suspended  particles  were  represented  by  streaks 
during  the  period  the  particles  were  within  the  illuminated 
sheet;  in  this  study  the  streaks  made  were  referred  to  as 
"pathlines. " 
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A.  LOW  POWER  OBSERVATIONS 


The  low  power  observation  was  made  at  a  3.2  W  power  level, 
corresponding  to  a  heat  flux  of  395  W/m2.  At  this  low  power, 
as  seen  in  Figure  7,  it  was  clear  that  meandering  of  the  plume 
was  taking  place.  In  the  top  photo  the  plume  flow  is  to  the 
right  and  seems  to  remain  steady  for  the  one  minute  exposure 
time.  The  bottom  photo,  taken  seconds  later  and  exposed  for 
two  minutes,  shows  the  plume  shifted  towards  the  left.  In 
this  photo  some  particles  are  still  flowing  to  the  right  and 
would  indicate  that  the  flow  to  the  left  was  not  quite  steady. 
The  second  plane  observations  at  this  power  level  were 
similar. 

Particles  in  both  planes  seemed  to  follow  distinct  paths 
except  when  the  plume  would  shift  while  meandering.  In  the 
plume  base  pathline  crossing  was  not  evident  and  would 
indicate  that  flow  at  this  low  power  was  laminar. 

Observations  at  a  slightly  higher  power  level  of  7.1  W 
(877  W/m2)  showed  a  distinct  flow  vertically  upward  with  no 
meandering  of  the  plume.  It  was  clear  from  the  visualization 
that  particles  were  still  following  pathlines  that  did  not 
cross.  Hence,  flow  was  assumed  to  be  laminar. 

B.  MEDIUM  POWER  OBSERVATIONS 

Observations  at  medium  power  levels  showed  characteristics 
of  departure  from  laminar  flow.  Photographs  were  taken  at 
power  levels  of  11.4  W  (1407  W/m2)  and  19.8  W  (2445  W/m2). 
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Figure  7.  Buoyant  Plume  at  Low  Power  (3.2  W)  (top)  60 
second  exposure  time  (bottom)  120  second  exposure  time. 
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Observations  were  also  made  for  a  pulsatile  input  power  using 
a  triangular  wave  pattern  with  an  average  power  of  19.8  W, 
0.1  Hz  frequency,  and  8  W  amplitude. 

Photographs  taken  for  the  nominally  steady  power  inputs 
clearly  showed  evidence  of  a  well  defined  buoyant  plume  with 
no  evidence  of  large  meandering  as  seen  at  lower  power  levels. 
As  seen  in  Figure  8  for  19.8  W  (2445  W/m2)  power  level, 
effects  of  time  dependence  are  observed  in  the  form  of 
crossing  of  the  particle  paths. 

Use  of  a  pulsatile  input  power  did  not  show  a  pronounced 
effect  on  the  buoyant  plume.  The  flow  pattern  for  a  19.8  W 
steady  input  power  was  not  significantly  different  from  that 
at  the  19.8  W  average,  pulsed  power  input.  Pathlines  for  both 
power  inputs,  as  shown  in  Figure  8,  are  quite  similar  within 
the  plume.  However,  outside  the  plume,  the  effect  of  the 
power  pulsation  at  this  low  frequency  appears  to  cause  a 
downward  flow  towards  the  heated  surface  prior  to  entrainment 
in  the  plume.  The  plume  width  is  also  seen  to  be  larger  in 
the  presence  of  pulsations.  This  is  obvious  in  the  bottom 
photograph.  Because  this  was  the  only  case  where  a  pulsatile 
power  was  used,  no  explanation  for  this  effect  is  possible. 

Also  evident  in  the  photographs  for  medium  power  levels 
was  a  slight  narrowing  of  the  plume  compared  to  the  lower 
input  power.  The  higher  power  level  causes  an  increase  in 
fluid  velocity  and  a  thinner  thermal  and  momentum  transport 
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region.  That,  coupled  with  the  appearance  of  crossing  of 
particle  paths  indicates  a  transition  to  turbulent  flow. 


C.  HIGH  POWER  OBSERVATIONS 

Observation  of  the  buoyant  flow  at  high  power  was  made  at 
48.9  W.  From  the  photograph  in  Figure  9,  a  further  narrowing 
of  the  plume  is  apparent.  Closer  to  the  surface,  pathlines 
appear  to  cross  much  more  than  in  the  photographs  at  lower 
powers.  Entrainment  of  particles  from  the  sides  along  the 
free  boundary  is  more  evident.  Flow  at  this  power  input 
appears  to  be  in  transition  or  turbulent  conditions  evidenced 
by  the  mixing  in  the  field  of  view  of  Figure  9. 


Figure  9.  Buoyant  Plume  at  High  Power  (48.9  w  -  60  second 
exposure  time) . 
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IV.  ANALYSIS  OF  PLUME  TEMPERATURE  DATA 

Plume  temperature  measurements  were  processed  in  two  ways. 
First,  the  temperature  difference  data  were  plotted  as  time 
series  graphs  as  shown  in  Figure  10.  The  x,y,  and  z 
coordinates  of  the  spatial  location  are  identified  at  the 
upper  right  corner.  A  total  of  54  measurements  were  taken  at 
each  power  level.  A  Fast  Fourier  Transform  (FFT)  was 
performed  to  better  identify  the  frequency  components  of  the 
temperature  signal  obtained  by  the  thermocouple  probe  over 
1024  data  points.  Data  were  initially  taken  at  a  15  Hz 
sampling  frequency  which  revealed  that  the  maximum  frequency 
observed  was  approximately  0.2  Hz.  Therefore,  the  sampling 
frequency  was  set  at  5  Hz,  in  excess  of  Nyquist  sampling 
criterion,  in  order  to  avoid  aliasing.  The  power  spectral 
density  (PSD) ,  a  measurement  of  the  energy  of  the  signal  at 
various  frequencies,  was  then  obtained.  The  frequencies  (in 
Hz)  defined  by  the  FFT  were  then  plotted  against  PSD  (a  non- 
dimensional  magnitude) ,  also  shown  in  Figure  10. 

Data  were  obtained  for  seven  constant  power  and  two 
periodic  power  inputs.  After  reviewing  the  plots  described 
above,  the  data  were  divided  into  three  categories:  low, 
medium,  and  high  power.  Next,  all  the  data  were  reviewed 
together  in  order  to  determine  group  trends.  Because  neither 
the  time  series  plots,  nor  the  FFT  plots  revealed  significant 
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differences  between  constant  power  inputs  and  the  periodic 
power  input  with  the  sane  mean,  periodic  power  inputs  will  not 
be  discussed. 


Time  (sec) 


Figure  10.  Example  of  Time  Series  and  FFT/PSD  Plots. 

A.  LOW  POWER  OBSERVATIONS 

Data  in  the  low  power  regime  were  obtained  at  0.78  w  (96 
W/m2)  and  3.2  W  (395  W/m2).  The  time  series  plots  for  the  0.78 
W  power  level  shown  in  Figure  11  reveal  low  frequency,  small 
amplitude  oscillations  close  to  the  surface,  and  are  nearly 
constant  at  all  other  elevations.  Seen  in  Figure  12,  the  FFT 
plots  at  this  power  level  show  only  a  few  peaks  with  small  PSD 


Power  Spectral  Density 


magnitudes  at  very  low  frequency.  Clearly,  buoyant  flow  at 
this  power  level  at  all  elevations  is  laminar. 

At  the  3.2  W  power  level  the  patterns  begin  to  change. 
Near  the  surface,  the  time  series  plot  shows  an  increase  in 
the  frequency  of  oscillation  with  larger  amplitude.  As  the 
elevations  increase,  the  dominant  frequency  of  oscillation 
increases,  and  the  amplitude  decreases.  See  Figures  13  and 
14.  The  FFT  plots  reveal  that  near  the  surface,  more 
frequencies  are  evident  with  much  higher  energy.  As  elevation 
increases,  more  frequencies  are  apparent,  but  with  lower 
energy.  At  the  lower  elevations,  the  buoyant  flow  is  clearly 
laminar.  But  as  elevation  increases,  it  is  believed  that  the 
flow  is  undergoing  transition  as  discussed  in  more  detail 
later . 

B.  MEDIUM  POWER  OBSERVATIONS 

Two  input  powers  were  considered  to  be  in  a  medium  regime 
because  both  the  time  series  and  FFT  plots  showed  trends  that 
were  noticeably  different  from  both  those  of  the  low  and  high 
power  regimes.  Medium  regime  data  were  obtained  at  12.2  W 
(1506  W/m2)  and  19.8  W  (2445  W/m2) .  For  both  power  levels, 
time  series  plots  showed  a  marked  increase  in  amplitude  and 
frequency  of  oscillation  of  delta  T  near  the  surface.  The  FFT 
plots  revealed  peaks  which  were  much  more  defined  and  over  a 
broader  range  of  frequencies.  The  energy  levels  at  the 
defined  frequencies  also  increased. 
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Figure  13.  Time  Series  along  Flume  Centerline  at  3.2  W 
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In  spite  of  the  similarities  mentioned  above,  differences 
between  the  12.2  H  and  19.8  W  power  levels  were  also  obvious. 
See  Figures  15  to  18.  The  time  series  plots  at  12.2  W  had 
smaller  amplitude  oscillations  as  compared  to  the  19.8  W  input 
power.  More  noticeable,  was  the  difference  in  the  FFT  plots. 
At  19.8  W,  many  more  frequencies  were  identified,  with  much 
higher  energy  levels.  It  is  believed  that  flow  near  the 
surface  at  12.2  W  is  laminar,  while  transitioning  at  all  other 
levels,  never  becoming  fully  turbulent.  The  broad  range  of 
frequencies  that  are  evident  in  the  FFT  plots  at  19.8  W  seems 
to  indicate  that  the  flow  at  the  two  elevations  near  the 
surface  is  in  transition,  as  shown  by  the  high-energy  low- 
frequency  peaks,  but  is  near  turbulent  at  higher  elevations 
where  more  lower  energy  frequencies  become  apparent.  As  was 
seen  in  the  (67.5,67.5)  location  plots  at  19.8  W,  a  broad 
range  of  frequencies  is  still  evident.  This  implies  that  the 
buoyant  plume  has  not  yet  narrowed,  an  indication  of  turbulent 
flow  evidenced  at  the  high  power  levels  in  flow  visualization. 

C.  HIGH  POWER  OBSERVATIONS 

Three  power  levels  were  considered  to  be  in  the  high  power 
regime:  27.8  W  (3432  W/m2)  ,  37.7  W  (4654  W/m2)  ,  and  48.9  W 
(6037  W/m2)  .  At  all  three  power  levels  the  buoyant  plumes 
have  narrowed.  This  is  clear  from  the  FFT  plots  at  the  off 
centerline  positions.  Of  the  three  power  levels,  at 
elevations  above  the  surface,  only  the  32  mm  elevation  at  37.7 
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Figure  15.  Time  Series  along  Flume  Centerline  at  12.2  W 
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Figure  16.  FFT/PSD  along  Plume  Centerline  at  12.2  w. 
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Figure  17.  Tin*  Series  along  Plume  Centerline  at  19.8  w. 
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Figure  18.  FFT/PSD  along  Pluae  Centerline  at  19.8  w. 
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W  resulted  in  a  noticeable  signal.  The  time  series  plots  at 
these  power  levels  show  a  nearly  constant  delta  T  over  the 
period  of  observation. 

At  the  plume  centerline  position,  time  series  plots  as 
well  as  the  FFT  plots  for  the  position  just  above  the  test 
surface,  show  a  trend  that  would  be  considered  to  be 
transition.  This  is  clear  from  the  FFT  plots  which  show  a 
variety  of  low  frequencies  at  very  high  energy.  At  the  higher 
elevations,  the  plots  would  indicate  that  the  flow  is 
transitioning  to  full  turbulent  flow  as  the  elevations 
increase.  At  all  power  levels  higher-energy  frequencies  are 
evident  over  the  broadest  spectrum.  See  Figures  19  to  24. 
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Tine  Series  along  Plume  Centerline  at  27.8  w. 
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Figure  21.  Tine  Series  along  Plune  centerline  at  37.7  V. 


Power  Spectral  Density 


o 

3000 


2000 

1000 


(45.45.32) 


i  i 


Lm 


3000 

2000 


JiL 


A 


1000 


(45.45,52) 


»■  - 


iA 


- -  S  ^ 


1000 


(45,45,82) 


M 


-£*L. 


2000 


1000  If! 


(45,45,122) 


^  --a 


VA- 


V— ftW 


(45,45,170) 


Ja_Vk ,  a A  a^v  _a  a 


01  07  03  fr  *  ,-  / 

o.l  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9 


Freq  (Hz) 

Figure  22.  FFT/PSD  along  Plume  Centerline  at  37.7  w. 
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Figure  23.  Time  Series  along  Plume  Centerline  at  48.9  w. 
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O.  COMBINED  ANALYSIS 

To  better  evaluate  the  transition  process  in  the  buoyant 
plume,  a  single  plot  was  developed  which  was  comprised  of  data 
from  all  the  FFT/PSD  graphs.  Data  for  the  various  x,y,z 
locations  in  the  plume  for  various  heat  fluxes  are 
represented.  A  spreadsheet  program  was  used  in  the 
development  of  the  plot,  with  physical  properties  of  the  water 
taken  at  a  film  temperature  defined  as  the  average  of  the  mean 
plate  and  ambient  temperatures,  and  neglecting  conduction  heat 
losses  across  the  bottom  of  the  test  surface.  An  elevation 
based  modified  Grashof  number: 


Gr 


* 

Z 


g§q"z* 

kv2 


was  plotted  against  a  non-dimensional  frequency: 


g$o" 

1/2 

k  J 

on  a  log-log  scale  as  shown  in  Figure  25.  Based  on  this  plot 
and  the  observations  made  from  FFT  data  and  flow 
visualizations,  conclusions  can  be  made  about  the  transition 
process . 

The  trends  of  the  data  points  represented  in  the  plot  can 
be  characterized  in  three  general  areas:  Gr*  <  500;  500  <  Gr* 
<  lxlO7;  and,  Gr*  >  lxlO7.  After  evaluating  all  the  graphs 
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mentioned  earlier  along  with  the  photographs  taken  during  flow 
visualization,  and  in  reviewing  this  final  plot,  it  is 
concluded  that  the  plot  should  be  divided  into  four  areas 
which  characterize  the  transition  process  based  on  power  and 
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Figure  25.  Non-dimensional  Frequency  versus  Grashof  Number. 

elevation  above  the  surface. 

l.  Laminar  Buoyant  Flow  Above  The  Surface 

The  first  region  to  consider  is  for  Gr*  <  500.  Data 
points  in  this  region  are  for  all  power  levels  at  an  elevation 
of  2  mm  above  the  surface.  As  reported  by  Bill  and  Gebhart 
[Ref.  12],  spectral  analysis  of  the  plume  in  the  laminar 
region  was  characterized  by  high-energy  frequencies  at  the  low 
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end  of  the  spectrum.  The  reason  for  the  appearance  of 
disturbances  is  attributed  to  swaying  and  oscillations  of  the 
plume  at  its  base.  This  seems  to  coincide  with  the  FFT/PSD 
plots  obtained. 

2.  Plume  Meandering 

For  the  region  500  <  Gr*  <  lxl 07  the  buoyant  flow 
above  the  surface  is  still  laminar  but  at  higher  elevations, 
the  plume  meanders.  Because  the  motion  is  a  slow  movement  of 
the  entire  plume  column,  no  oscillations  are  detected  in  the 
spectrum  analysis.  The  meandering  motion  is  only  evident  for 
low  power  levels. 

3.  Plume  Transition  Region 

Transition  of  the  Plume  is  believed  to  take  place  for 
lxio7  <  Gr*  <  6xl09.  The  beginning  point  for  transition  was 
concluded  to  be  somewhat  higher  than  the  values  plotted  for 
the  3.2  W  power  level  at  the  32  mm  elevation.  From  the 
photographs  obtained  at  this  power  level,  plume  meandering  was 
observed.  The  end  of  transition  was  based  on  the  high  end 
data  points  for  12.2  W.  Clearly  from  the  FFT/PSD  plots  for 
the  centerline  position  at  this  power  level,  the  frequencies 
represented  do  not  cover  a  wide  spectrum  characterizing  fully 
turbulent  conditions.  Also,  the  time  series  plots  do  not 
represent  a  highly  oscillating  delta  T. 
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4.  Turbulent  Buoyant  Flow 

For  Gr*  £  6xl09  flow  is  concluded  to  be  turbulent. 
From  FFT/PSD  plots  of  the  data  points  represented  by  the  power 
levels  in  this  region,  frequencies  of  substantial  energy  are 
apparent  across  a  large  spectrum.  Many  data  points  for  the 
higher  power  levels  were  not  utilized  in  Figure  25  because  so 
many  distinct  frequencies  were  apparent  that  the  flow  was 
correctly  assumed  turbulent.  Hence,  the  lack  of  data  points 
in  the  turbulent  region. 
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V.  TEST  SURFACE  TEMPERATURE  ANALYSIS 
A.  CONSTANT  POWER 

Figure  26  shows  the  heated  surface  centerline  temperature 
pattern  for  six  of  the  seven  power  levels  tested.  Assuming 
that  conduction  losses  in  the  surface  are  negligible,  the 
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Figure  26.  Heater  Surface  Temperatures  at  constant  Power. 

power  supplied  to  the  heater  is  equal  to  the  heat  transferred 
to  the  water  by  convection.  Non-dimensional  heat  transfer 
data  were  based  on  the  net  heat  input  rates . 
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The  Nusselt  (Nu)  and  Rayleigh  (Ra)  numbers  were  calculated 


as: 


*r  _  _ Gconv^c _ 

“  ■  kAiT^-T^) 


Ra 


3$  Qconv^c 

kAv2 


Where  Tplilc  is  the  average  plate  temperature,  Lc  is  a 
characteristic  length  taken  as  the  ratio  of  the  heater  surface 
area  to  perimeter . [Ref .  8] 

Having  taken  all  fluid  properties  at  a  film  temperature 
defined  as  the  average  of  the  plate  and  ambient  fluid 
temperatures,  a  plot  of  the  Nusselt  number  as  a  function  of 
flux  based  Rayleigh  number  was  made.  Figure  27  shows  a  linear 
best  fit  of  the  plotted  data  which  resulted  in  the  following 
correlation: 

Nu  =  0.135  Ra0227 

valid  for  7.5xl05  <  Ra  <  8.3xl07.  This  shows  a  higher  exponent 
of  Ra  than  the  correlation  determined  by  Hickey: 

Nu  =  2.22  Ra0  065 

which  was  valid  for  1.2xl06  <  Ra  <  9.1xl07.  This  could  be 
caused  by  a  greater  thermal  resistance  between  the  heater  and 
measuring  thermocouples  in  Hickey's  data.  Evaluation  of  a 
temperature  based  Rayleigh  number  as  defined  by  Hickey 
permitted  comparison  to  correlations  found  by  Lloyd  and  Moran 


48 


Figure  27.  Average  Musselt  Number  versus  Flux  Based 
Rayleigh  Number. 

[Ref.  4].  A  linear  best  fit  of  the  data  in  this  plot  shown  in 
Figure  28  resulted  in  the  following  correlation: 

Nu  =  0.08  Ra0293 

valid  for  2.6x10s  <  Ra  <  9.8xl06.  As  seen  in  the  figure,  the 
slope  of  the  present  data  agrees  with  that  of  Lloyd  and  Moran 
but  the  actuhl  values  are  considerably  lower.  Hickey's  data 
are  closer  but  display  a  lower  slope,  perhaps  due  to  a  higher 
thermal  resistance  between  the  plate  and  the  measuring 
thermocouples  as  noted  earlier.  As  pointed  out  by  Hickey, 
Lloyd  and  Moran  developed  their  correlations  based  on  a  fully 
heated  test  surface  while  the  present  study  has  a  large 
unheated  low  thermal  conductivity  substrate;  which  impedes 
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Figure  28.  Average  Nusselt  Number  versus  Temperature  Based 
Rayleigh  Number. 

heat  transfer  near  the  edges.  Also,  the  fluid  used  by  Lloyd 
and  Moran  was  air  while  the  present  study  was  based  on  water 
and  differences  in  the  fluid  Prandtl  number  may  effect  heat 
transfer.  Lastly,  Lloyd  and  Moran  developed  separate 
correlations  for  laminar  and  turbulent  flow.  The  correlations 
developed  by  Hickey  and  that  of  the  present  study  are  for  the 
entire  range  of  Rayleigh  numbers  examined.  The  correlation 
developed  by  Lloyd  and  Moran  for  laminar  flow  had  a  slope  of 
0.255  and  was  valid  for  2.6xl04  <  Ra  <  8xl06  which  is  in  fair 
agreement  with  the  present  study. 
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B.  PERIODIC  INPUT  POWER 

As  was  described  in  the  previous  chapter,  the  present 
study  showed  no  apparent  benefit  in  using  a  periodic  input 
power.  Hickey  [Ref  8.]  noted  that  for  certain  types  of  input 
powers  at  specific  frequencies,  enhancement  in  the  heat 
transfer  from  the  surface  could  be  obtained.  Several  attempts 
were  made  to  repeat  results  obtained  by  Hickey,  but  no 
significant  heat  transfer  enhancement  was  noted.  From  flow 
visualization  it  was  apparent  that  the  flow  was  altered  by  the 
pulsatile  power  input.  However,  too  little  data  was  obtained 
in  the  present  study  for  pulsatile  input  powers  to  verify  the 
heat  transfer  enhancement  findings  made  by  Hickey. 
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VI .  CONCLUSIONS 


Use  of  a  traversing  thermocouple  probe  coupled  with  flow 
visualization  proved  invaluable  in  conceptualizing  the  flow 
patterns  within  a  buoyant  plume.  The  temperature  measurements 
obtained  by  the  probe  were  significant  in  determining  the 
extent  of  the  plume  over  the  test  surface.  Flow  visualization 
was  vital  in  verifying  the  meandering  of  the  plume.  Both 
programs  were  equally  important  in  understanding  the 
transition  process  for  a  free  boundary  flow  above  a  flat 
horizontal  heated  surface  surrounded  by  an  unheated  area. 
This  can  contribute  significantly  to  studies  in  liquid 
immersion  cooling,  particularly  in  the  validation  of 
computational  codes. 

Use  of  a  pulsatile  input  power  provided  no  effect  in 
obtaining  heat  transfer  enhancement.  The  FFT/PSD  plots  for 
the  pulsatile  input  powers  revealed  that  the  frequency  of  the 
input  power  was  buried  in  the  temperature  signal  obtained.  It 
is  believed  that  the  0.1  Hz  input  power  frequency  was  too  low 
to  affect  the  plume. 

The  reconstructed  test  surface  seemed  to  work  as 
effectively  as  that  built  by  Hickey  [Ref. 8].  However,  the 
heat  transfer  correlations  developed  from  the  data  obtained 
from  the  present  surface  seemed  to  be  in  better  agreement  with 
the  trends  of  the  previous  studies. 
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VII.  RECOMMENDATIONS 

While  the  present  study  was  a  continuation  of  that  begun 
by  Hickey  [Ref.  8],  further  research  would  make  the  study 
complete.  Therefore,  it  is  recommended  that  the  following 
aspects  be  studied  further: 

•  Use  the  data  and  results  presented  in  this  study  to 
develop  analytical  models  to  predict  the  transition 
process  of  buoyant  plumes. 

•  The  photographs  obtained  for  through  flow  visualization 
were  very  limited.  Flow  visualization  should  be  repeated 
to  obtain  photographs  chat  clearly  show  the  flow  from  the 
test  surface  as  well  as  at  a  higher  level  in  the  buoyant 
plume . 

•  Determine  the  velocity  field  in  and  around  the  plume  using 
a  laser  doppler  velocimeter  to  better  understand  the  flow 
patterns . 

•  Develop  theoretical  models  which  can  utilize  the  results 
obtained  in  order  to  better  determine  the  heat  transfer 
effects  on  the  heater  surface. 
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APPENDIX  A 


SAMPLE  CALCULATIONS 

The  following  calculations  were  made  for  a  heat  flux  of 
3432  W/m2,  f  =  0.425  Hz,  AT  =  17.42°C,  and  at  the  elevation 
z  =  0.052  m. 

1.  Characteristic  Dimensions. 

Perimeter  (P)  =  4  x  0.09  m 

=  0.36  m 

Heater  Surface  Area  (A,)  =  (0.09  m)2 

=  0.0081  m2 

Characteristic  Length  (Lc)  =  A,/P 

=  0.0225  m 

2.  Convective  Heat  Flux. 

Input  Power  (Q^)  =  27.8  W 
Heat  Flux  (q")  =  27.8/0.0081 
=  3432  W/m2 

3.  Fluid  Properties  (Ref.  17]. 

In  the  plume: 

0  =  263.3X10-6  K'1 
v  =  0.885X10-6  m2/s 
k  =  0.611  W/m-K 
Pr  =  5.985 
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On  the  surface: 

/3  =  334.6x10-*  K'* 
v  =  0.749x10-*  m2/s 
k  =  0.623  W/rn-K 
Pr  =  5.00 

4.  Mon-dimensional  Frequency. 

f*  =  0 . 425/ [9.81  *  263.3E-6  *  3432/0 . 611] * 
=  0.111 

5.  Modified  Grashof  Number. 

G  •  =  9.81  •  263 . 3E-6  •  3432  •  0 . 052* 

*  0.611  •  (0.885E-6)2 

Gl*z  -  1.36 X106 

6.  Musselt  Number. 

Nu  =  [3432  •  0 . 0225 ] / [0 . 623  *  17.42] 

=  7.12 

7 .  Rayleigh  Numbers . 

Flux  based: 

Ra  =  9 -81  ‘  334 ■ 6E-6  •  3432  •  0.022545  Q 
0.623  *  (0.749F-6)2 

Ra  =  4.13 X101 
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Temperature  based 


Ra 


9.81  •  334.6 E-6  ■  17.42  • 
(0 . 749JJ-6 )  2 


Ra  =  5 . 8x1 06 


0,02_25ZS0 
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APPENDIX  B 


UNCERTAINTY  ANALYSIS 


The  accuracy  of  the  data  in  this  study  was  determined  by 
performing  an  uncertainty  analysis.  For  a  function  comprising 
a  number  of  independent  measurements  F  =  F(XA,  XB,  Xc)  ,  the 
uncertainty  of  F  was  calculated  as: 


6F  = 


If  the  function  F  =  C  XA*  XBb  X^ ,  then  the  uncertainty  is 
determined  by: 


bF 

F 


List  of  properties  [Ref.  17]: 
g  =  9.81  m/s2 

Tflto  =  298  K:  Tflta  =  306.8  K: 

/8  =  263 . 3X10"6  K'1  0  =  334.6x10^  K1 

k  =  0.611  W/ro-K  k  =  0.623  W/m  K 

v  =  0.885X10'6  m2/s  v  =  0 . 749x10"*  m2/s 

Pr  =  5.985  Pr  =  5.000 
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1.  Frequency  Uncertainty 


To  convert  the  time  series  to  the  frequency  domain,  a  Fast 
Fourier  Transform  was  performed  over  1024  data  points.  The 
highest  measurable  frequency  was  determined  as: 

f  ,  _ i _ 

max  2  Period  of  Time  Record 


Therfore,  f„„  =  1.97  Hz  and  fif=0.002  Hz.  Hence, 


Ml 

f 


|l/2 


for: 


bq"  = 


and  with: 

Q  =  27.8  W 
A  =  0.0081  m2 

then, 

q"  =  3432  W/m2 


S Q  =  0.0001  W 
Sk  =  3.6X10'7  m2 


S q"  =  42.4  W/m2 
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With, 

f  =  0.425  Hz  <Sf  =  0.002  Hz 

The  uncertainty  of  S f*/f*  =  0.0078  or  0.78% 

2.  Modified  Grashof  Number  Uncertainty. 


Using  the  values  for  q"  and  £q"  from  section  1  and  with  z  = 
0.052  m  and  Sz  -  0.0001  m,  then  £Grt*/Grz*  -  0.0145  or  1.45%. 

3.  Nusselt  Number  Uncertainty. 


AT 

SLC 


TpUlr  “  ^mk 

2.78X10"* 


II 

S’ 

•o 

&g"f  + 

/  6At\2 

1/2 

Nu  [\ 

q"  ) 

{  Lc  I 

\  Ar) 

=  17.42° 

C  and 

SAT  -  0. 

1°C, 

Lc 

m,  then 

5Nu/Nu 

=  0.0184 

or  1 . 

84% 

=  0.0225  m  and 
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4.  Rayleigh  Number  Uncertainty 


Ra 


9fiQ"L*c 

kv2 


Pi 


6  Ra 
Ra 


Then  £Ra/Ra  =  0.051  or  5.1%. 
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